The Role of Subsurface Flows in Solar Surface Convection: Modeling the Spectrum of Supergranular and Larger Scale Flows

We model the solar horizontal velocity power spectrum at scales larger than granulation using a two-component approximation to the mass continuity equation. The model takes four times the density scale height as the integral (driving) scale of the vertical motions at each depth. Scales larger than this decay with height from the deeper layers. Those smaller are assumed to follow a Kolomogorov turbulent cascade, with the total power in the vertical convective motions matching that required to transport the solar luminosity in a mixing length formulation. These model components are validated using large scale radiative hydrodynamic simulations. We reach two primary conclusions: 1. The model predicts significantly more power at low wavenumbers than is observed in the solar photospheric horizontal velocity spectrum. 2. Ionization plays a minor role in shaping the observed solar velocity spectrum by reducing convective amplitudes in the regions of partial helium ionization. The excess low wavenumber power is also seen in the fully nonlinear three-dimensional radiative hydrodynamic simulations employing a realistic equation of state. This adds to other recent evidence suggesting that the amplitudes of large scale convective motions in the Sun are significantly lower than expected. Employing the same feature tracking algorithm used with observational data on the simulation output, we show that the observed low wavenumber power can be reproduced in hydrodynamic models if the amplitudes of large scale modes in the deep layers are artificially reduced. Since the large scale modes have reduced amplitudes, modes on the scale of supergranulation and smaller remain important to convective heat flux even in the deep layers, suggesting that small scale convective correlations are maintained through the bulk of the solar convection zone.Comment: 36 pages, 6 figure

Suppression of Kelvon-induced decay of quantized vortices in oblate Bose-Einstein Condensates

We study the Kelvin mode excitations on a vortex line in a three-dimensional trapped Bose-Einstein condensate at finite temperature. Our stochastic Gross-Pitaevskii simulations show that the activation of these modes can be suppressed by tightening the confinement along the direction of the vortex line, leading to a strong suppression in the vortex decay rate as the system enters a regime of two-dimensional vortex dynamics. As the system approaches the condensation transition temperature we find that the vortex decay rate is strongly sensitive to dimensionality and temperature, observing a large enhancement for quasi-two-dimensional traps. Three-dimensional simulations of the recent vortex dipole decay experiment of Neely et al. [Phys. Rev. Lett. 104, 160401 (2010)] confirm two-dimensional vortex dynamics, and predict a dipole lifetime consistent with experimental observations and suppression of Kelvon-induced vortex decay in highly oblate condensates.Comment: 8 pages, 8 figure

The Water holes at Ijara

Volume: XXI

Stratified shear flow instabilities at large Richardson numbers

Numerical simulations of stratified shear flow instabilities are performed in two dimensions in the Boussinesq limit. The density variation length scale is chosen to be four times smaller than the velocity variation length scale so that Holmboe or Kelvin-Helmholtz unstable modes are present depending on the choice of the global Richardson number Ri. Three different values of Ri were examined Ri =0.2, 2, 20. The flows for the three examined values are all unstable due to different modes namely: the Kelvin-Helmholtz mode for Ri=0.2, the first Holmboe mode for Ri=2, and the second Holmboe mode for Ri=20 that has been discovered recently and it is the first time that it is examined in the non-linear stage. It is found that the amplitude of the velocity perturbation of the second Holmboe mode at the non-linear stage is smaller but comparable to first Holmboe mode. The increase of the potential energy however due to the second Holmboe modes is greater than that of the first mode. The Kelvin-Helmholtz mode is larger by two orders of magnitude in kinetic energy than the Holmboe modes and about ten times larger in potential energy than the Holmboe modes. The results in this paper suggest that although mixing is suppressed at large Richardson numbers it is not negligible, and turbulent mixing processes in strongly stratified environments can not be excluded.Comment: Submitted to Physics of Fluid

Stochastic time-dependent current-density functional theory: a functional theory of open quantum systems

The dynamics of a many-body system coupled to an external environment represents a fundamentally important problem. To this class of open quantum systems pertains the study of energy transport and dissipation, dephasing, quantum measurement and quantum information theory, phase transitions driven by dissipative effects, etc. Here, we discuss in detail an extension of time-dependent current-density-functional theory (TDCDFT), we named stochastic TDCDFT [Phys. Rev. Lett. {\bf 98}, 226403 (2007)], that allows the description of such problems from a microscopic point of view. We discuss the assumptions of the theory, its relation to a density matrix formalism, and the limitations of the latter in the present context. In addition, we describe a numerically convenient way to solve the corresponding equations of motion, and apply this theory to the dynamics of a 1D gas of excited bosons confined in a harmonic potential and in contact with an external bath.Comment: 17 pages, 7 figures, RevTex4; few typos corrected, a figure modifie

Suppressing the Rayleigh-Taylor instability with a rotating magnetic field

The Rayleigh-Taylor instability of a magnetic fluid superimposed on a non-magnetic liquid of lower density may be suppressed with the help of a spatially homogeneous magnetic field rotating in the plane of the undisturbed interface. Starting from the complete set of Navier-Stokes equations for both liquids a Floquet analysis is performed which consistently takes into account the viscosities of the fluids. Using experimentally relevant values of the parameters we suggest to use this stabilization mechanism to provide controlled initial conditions for an experimental investigation of the Rayleigh-Taylor instability

Periodic magnetorotational dynamo action as a prototype of nonlinear magnetic field generation in shear flows

The nature of dynamo action in shear flows prone to magnetohydrodynamic instabilities is investigated using the magnetorotational dynamo in Keplerian shear flow as a prototype problem. Using direct numerical simulations and Newton's method, we compute an exact time-periodic magnetorotational dynamo solution to the three-dimensional dissipative incompressible magnetohydrodynamic equations with rotation and shear. We discuss the physical mechanism behind the cycle and show that it results from a combination of linear and nonlinear interactions between a large-scale axisymmetric toroidal magnetic field and non-axisymmetric perturbations amplified by the magnetorotational instability. We demonstrate that this large scale dynamo mechanism is overall intrinsically nonlinear and not reducible to the standard mean-field dynamo formalism. Our results therefore provide clear evidence for a generic nonlinear generation mechanism of time-dependent coherent large-scale magnetic fields in shear flows and call for new theoretical dynamo models. These findings may offer important clues to understand the transitional and statistical properties of subcritical magnetorotational turbulence.Comment: 10 pages, 6 figures, accepted for publication in Physical Review

Cavitation and bubble collapse in hot asymmetric nuclear matter

The dynamics of embryonic bubbles in overheated, viscous and non-Markovian nuclear matter is studied. It is shown that the memory and the Fermi surface distortions significantly affect the hinderance of bubble collapse and determine a characteristic oscillations of the bubble radius. These oscillations occur due to the additional elastic force induced by the memory integral.Comment: Revtex file (10 pages) and 3 figure

Polarizability of conducting sphere-doublets using series of images.

The classical electrostatic problem of two nonintersecting conducting spheres in a uniform incident electric field is considered. Starting from the basic Kelvin’s image principle, the two spheres are replaced with equivalent series of image sources, from which the polarizability is calculated. Explicit expressions for the axial and transversal components of the polarizability dyadic are found by solving the recurrence equations. Efficient numerical evaluation of the different series is also discussed.Peer reviewe

The shape and erosion of pebbles

The shapes of flat pebbles may be characterized in terms of the statistical distribution of curvatures measured along their contours. We illustrate this new method for clay pebbles eroded in a controlled laboratory apparatus, and also for naturally-occurring rip-up clasts formed and eroded in the Mont St.-Michel bay. We find that the curvature distribution allows finer discrimination than traditional measures of aspect ratios. Furthermore, it connects to the microscopic action of erosion processes that are typically faster at protruding regions of high curvature. We discuss in detail how the curvature may be reliable deduced from digital photographs.Comment: 10 pages, 11 figure